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Abstract 

The calculations in Thomas-Fermi approximation show that in a 
gravitational field each cell of ultra dense matter inside celestial bodies 
obtains a very small positive electric charge. A celestial body is elec- 
trically neutral as a whole, because the negative electric charge exists 
at its surface. The positive volume charge is very small, on the order 
of magnitude it equals to 10 _18 e per atom only. But it is sufficient to 
explain the occurrence of magnetic fields of the celestial bodies and 
the existence of a discrete spectrum of steady-state values of masses 
of planets, stars, and pulsars. 
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1 Introduction 



According to the conventional point of view, gravity does not induce any 
electric polarization in the interior of celestial bodies and electric forces are 
never considered in the balance of matter of celestial bodies. Moreover, it 
is generally assumed that the electric interaction plays practically no role in 
astrophysics. It is a consequence of the comprehension that the apprecia- 
ble electric polarization cannot arise in metals and other nonsegneto- and 
nonpiro-electric materials. It is entirely correct for all substances under ac- 
tion of small pressure. But, thus, one can disregard the fact that ultrahigh 
pressure transmutes all substances into plasma state and radically changes 
the properties of substance. In ultradense plasma, there is a different addi- 
tional mechanism of the gravity-induced electric polarization. 

In a large celestial body, consisting of ultradense plasma, this gravity- 
induced electric polarization (GIEP) can be rather great and can play a 
determining role in the formation of a number of features of the structure of 
a celestial body and its properties. 

First of all, it concerns the following three problems, the statement and 
the solution of which change drastically: 

- the distribution of pressure and density of matter inside a celestial body; 

- the generation of a magnetic field by celestial bodies; 

- the formation of a spectrum of steady-state values of masses of celestial 
bodies. 

As a consequence, these features of the structure can influence the evolu- 
tion of stars. 



2 The gravity-induced electric polarization in 
conducting matter 

The action of gravity on metals has often been a topic of discussion before 
[|]]-||. The basic result of these researches is reduced to the statement that 
inside a metal gravity induces an electric field with an intensity 



B=S55, (i) 



where m,- is the mass of an ion, 
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g is gravity acceleration, 
e is the electron charge. 

This field is so small that it is not possible to measure it experimentally. 
It is a direct consequence of the presence of an ion lattice in a metal. This 
lattice is deformed by gravity and then the electron gas adapts its density to 
this deformation. The resulting field becomes very small. 

Under superhigh pressure, all substances transform into ultradense matter 
usually named nuclear-electron plasma 0. It occurs when external pressure 
enhances the density of matter several times [0, ||. Such values of pressure 
exist inside celestial bodies. 

In nuclear-electron plasma the electrons form the degenerated Fermi gas. 
At the same time, the positively charged ions form inside plasma a dense 
packing lattice [pf , [fL0|l . As usually accepted, this lattice may be replaced 
by a lattice of spherical cells of the same volume. The radius r s of such a 
spherical cell in plasma of the mass density 7 is given by 

^j-m- 1 -*, (2) 

3 \rriiJ n 

where Z is the charge of the nucleus, m; = Am p is the mass of the nucleus, 
A is the atomic number of the nucleus, m p is the mass of a proton, and n is 
the electron number density 

n = — (3) 
47rr^ 

The equilibrium condition in matter is described by the constancy of 
its electrochemical potential J7|. In plasma, the direct interaction between 
nuclei is absent, therefore the equilibrium in a nuclear subsystem of plasma 
(at T = 0) looks like 

Hi = rriiip + Zap = const. (4) 

Here <p is the potential of an electric field and ip is the potential of a 
gravitational field. 

The direct action of gravitation on electrons can be neglected. Therefore, 
the equilibrium condition in the electron gas is 

He = 7^- - (e - 5q)(p = const, (5) 
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where m e is the mass of an electron and pf is the Fermi momentum. 

By introducing the charge Sq, we take into account that the charge of the 
electron cloud inside a cell can differ from Ze. A small number of electrons 
can stay at the surface of a plasma body where the electric potential is 
absent. It results that the charge in a cell, subjected to the action of the 
electric potential, is effectively decreased on a small value Sq. If the radius 
of a star R is approximately 10 10 cm, one can expect that this mechanism 
gives on the order of magnitude ^ ~ ~ 10~ 18 . 

The electric polarization in plasma is a result of changing in density of 
both nuclear and electron gas subsystems. The electrostatic potential of the 
arising field is determined by the Gauss' law 



1 d 

r 2 dr 



o d 




r —(p 
dr 


= -47T 





ZeS{r) 



en 



(6) 



where the position of nuclei is described by the function 5(r). 
According to the Thomas - Fermi method, n is approximated by 



n 



8vr 
3^3 



Pf- 



(7) 



With this substitution, Eq.(|6|) is converted into a nonlinear differential 
equation for ip, which for r > is given by 



1 d ( 2 d 
r z dr \ dr 



Air 



in 



3h 3 



[2m e (/i e + (e-<5gV)] 3/2 



It can be simplified by introducing the following variables [|TT 

/i e + (e — Sq)f = Ze 2 — 

r 

and r = ax, 
where 

_ r 9?r 2 1 1/3 „ 
a - 1T28ZJ °0 

with an = -^-2 = Bohr radius. 
With the account of Eq. (H) 



(8) 



(9) 



Ze — = const 



niiip 



6q<p. 



(10) 



Then Eq.® gives 
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d 2 u w 3//2 
dx 2 x 1 / 2 



(11) 



In terms of u and x, the electron density within a cell is given by [[I I 



57T 



32Z 2 fu\ 3 / 2 



(12) 



Under the influence of gravity, the charge of the electron gas in a cell 
becomes equal to 



Q e = 4ne / n(r)r dr = — — 
Jo in 6 

Using Eq.(|ll|), we obtain 



2m, 



Ze 2 



3/2 



Ana 6 / x dx 
Jo 



u 

IX 



3/2 



(13) 



Q, = Z< j xdx^-^ = Ze [ dx 



o dx i dx 



du 
x— u 



Ze 



x. 



du 
dx 



-u(x s )+u(Q) 

(14) 

At r — > the electric potential is due to the nucleus alone ip(r) — > — . It 
means that w(0) — * 1 and each cell of plasma obtains a small charge 



Sq = Ze 



du 
dx 



- u(x. 



Zex s 



d / u 
doc \ oc 



For a cell placed in the point R inside a star 

dR 



5q = Zer\ 



d (u 
dR\r 



dr. 



(15) 



(16) 



Considering that the gravity acceleration g = — ^| and the electric field 



intensity E = - J| 



dr* ri 



5q 



dTl e 

This equation has the following solution 

dr 



dR 







(17) 



(18) 



5 



and 



^g + 5qE = 0. (19) 

In plasma, the equilibrium value of the electric field on nuclei according 
to Eq.([|) is determined by Eq.(|I]) as well as in a metal. But there is one more 
additional effect in plasma. Simultaneously with the supporting of nuclei in 
equilibrium, each cell obtains an extremely small positive electric charge. 

As divg = —AirGnmi and divE = 4im5q, the gravity-induced electric 
charge in a cell 

Sq = VG^ ~ l(T 18 e, (20) 
where G is the gravity constant. 

However, because the sizes of bodies may be very large, the electric field 
intensity may be very large as well 

E=JL. (21) 



In accordance with Eqs.(|IB|,19), the action of gravity on matter is com- 
pensated by the electric force and the gradient of pressure is absent. 

Thus, a celestial body is electrically neutral as a whole, because the pos- 
itive volume charge is concentrated inside the charged core and the negative 
electric charge exists on its surface and so one can infer gravity-induced elec- 
tric polarization of a body. 



3 Pressure distribution inside a celestial body. 

As at the surface of a celestial body pressure is absent, near this surface 
there is always a stratum where plasma and polarization are absent. For the 
large stars, the size of this stratum is insignificant. But for a small planet it 
can comprise a substantial part of a planet, and thus, only a small relatively 
internal region will be polarized. At the surface of this core, the electric field 
intensity falls to zero. The jump in the electric field intensity is accompanied 
on the surface of the core by the pressure jump Ap(i2jv)(l3"E31)- ^he 



important astrophysical consequence of the GIEP effect is the redistribution 
of the matter density inside a celestial body. In a celestial body, consisting 
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of matter with an atomic structure, density and pressure grow monotonously 
with depth. In a celestial body, consisting of electron-nuclear plasma, the 
GIEP effect results in the fact that the pressure gradient inside the polarized 
core is absent and the matter density is constant. Pressure affecting the 
matter inside this body is equal to the pressure jump at the surface of the 
core 

P = MRn) = = ^G 7 2 i4, (22) 

where 7 is the matter density in the core and Rat is the radius of the core. 

One can say that this pressure jump is due to the existence of the po- 
larization jump or, which is the same, the existence of the bonded surface 
charge, which is formed by electron pushed out from the core and which 
makes the total charge of the celestial body equal to zero. 



4 Earth's structure. 

It is important, that the GIEP effect gives the possibility to construct the 
intrinsically self-consistent theory of the Earth Although it is rather a 
solution of a geophysical problem than an astrophysical effect. 

Earlier models of the Earth assumed the existence of the monotonous 
dependence of pressure inside the planet. The division of the Earth into the 
core and the mantle was explained by the fact that at the creation of the 
Earth, on its share a certain amount of iron (and other heavy metals) and 
also a necessary amount of stone were given out. The core consists of metals 
and the mantle consists of stone. In these models, it was necessary to fit the 
parameters to get the densities of core and mantle and their sizes. It is not 
necessary to introduce any free parameters into the Earth theory based on 
the GIEP effect. Assuming that the Earth consists of homogeneous matter, 
the division on core and mantle is explained by the existence of the pressure 
jump on the surface of the core Eq.(^). The basic results of this theory are 
reduced to the calculation of the following five values: 

a) the radius of the Earth's core; 

b) the density of core matter; 

c) the density jump on the core-mantle boundary; 

d) the mass related to one electron of the Fermi gas in the core; 
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e) the electric polarization of the core. 

To express it in appropriate equations, one should substitute the following 
four parameters (the gravitational constant G is known): 

a) the mass of the Earth; 

b) the radius of the Earth; 

c) the matter density on the surface of the Earth; 

d) the bulk module of matter at the surface of the Earth. 

Thus, other parameters can be obtained, for example, the pressure dis- 
tribution inside the Earth. The basic results of this theory are shown in 
Fig.l. 

In addition, from the obtained data it is possible to calculate the angular 
momentum of the Earth. This calculation gives the value of 0.339MR 2 . It 
is in agreement with the measured value of 0.331MR 2 within several percent 
of the accuracy. 

It is possible to calculate the magnetic moment of the Earth. 

Apparently, using the appropriate data of other planets (the mass, the 
size, and the properties of matter at the surface), it is possible to construct 
models of these planets. It can be made, if these planets have electrically 
polarized cores and corresponding magnetic fields. 



5 The gyromagnetic ratio of a celestial body 

Another astrophysical consequence of the GIEP effect is coupled by the ro- 
tation of celestial bodies about their axes. A celestial body is electroneutral 
as a whole. The positive volume charge is concentrated inside the core and 
the negative charge is located at the surface of the core. When rotating, 
they move on different radii. As a result, all celestial bodies, when the GIEP 
effect is present, obtain magnetic moments 

2 4ir 

V = ^P QR n- ( 23 ) 

If the size of the body is sufficiently large, the core radius Rn does not 
differ significantly from its external radius R. For this celestial body, the 
angular momentum of the core coincides by the order of magnitude with the 
angular momentum of the body as a whole 
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Figure 1: The radial dependence of pressure and the matter density inside the 
Earth. The solid line is the calculated dependence of the matter density; the 
dashed line is the density of the Earth obtained by measuring the propagation 
velocity of seismic waves. The dash-dotted line is the calculated dependence 
of pressure inside the Earth over bulk module B=1.3 • 10 12 dyn/cm 2 . 
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L = -MttR 2 
5 



(24) 



where M = ^-jR 3 is the mass of a celestial body and Q is the velocity of 
rotation. 

Finally, the gyromagnetic ratios for these bodies should be close to the 
universal value 



L 



G 1 ' 2 
3c 



(25) 



The values of /x(L) for all celestial bodies (for which they are known 
today) are shown in Fig.2. The data for planets are taken from [14|, the data 
for stars are taken from 
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16|. 



and for pulsars - from 
As can be seen from the figure with the logarithmic accuracy, all celestial 
bodies - stars, planets, and pulsars - really have the gyromagnetic ratio close 
to the universal value ^r-. Only the data for the Moon fall out, because its 



3c 

size and inner pressure are too small to create an electrically polarized core. 
The estimation of the magnetic moment of the Earth within the frame of the 
theory mentioned above [13] gives fi ~ 4 • 10 25 Gs • cm 3 . It is almost precisely 
one half from the observed value of 8.05 • 10 25 Gs • cm 3 . For some planets, 
the values of magnetic moments are in a good agreement with Eq.(^) but 
they have an opposite sign. Apparently, it means that the hydrodynamic 
mechanism also plays a certain role. 

For the majority of pulsars, there are estimations of magnetic fields [H5 
obtained using a number of model assumptions [16]. It is impossible to 



consider these data as the data of measurements, but nevertheless, they also 
agree in certain way with Eq.(p5|),(Fig.3) 



6 The masses of celestial bodies. 

The important astrophysical outcome of the GIEP effect is a discrete distri- 
bution of masses of celestial bodies. This spectrum is a result of the fact that 
electron-nuclear plasmas can exist in various states. 

The equation of state of matter subjected to high pressure is usually 
described as a polytrope |7|]: 
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Figure 2: The observed values of the magnetic moments of celestial bodies 
vs. their angular momenta. On the ordinate, the logarithm of the magnetic 
moment over Gs ■ cm 3 is plotted; on the abscissa the logarithm of the angular 
momentum over erg ■ s is shown. The solid line illustrates Eq.(^). The dash- 
dotted line is the fitting of the observed values. 
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p = C-7 1+ *, (26) 

where C is the dimensional constant, 
k is the polytropy. 



6.1 Nonrelativistic electron-nuclear plasma. 

At relatively small pressure, substances are transmuted into nonrelativistic 
electron-nuclear (or electron-ion) plasma. It is peculiar to conditions existing 
inside cores of planets. According to @, the state equation of the nonrela- 
tivistic electron-nuclear plasma (characterized by the polytropy k=3/2) is 

_ (3^) 2/3 

Pi3/2) ~ 5m e (f3-m p )^ (2?) 

where (3 ■ m p is the mass of matter related to one electron of the Fermi 
gas system and m p is the proton mass. 

If the pressure inside a celestial body is formed by the GIEP effect and is 
determined by Eq . fl22|) , than from Eq. (|27|) for the nonrelativistic Fermi gas 
of electrons, we obtain the steady-state value of mass for a core of planet 

/ h 2 \ 3/2 y/a 

M (3 /2) = C (3 /2) • tt— — • , (28) 



Gm e m p I (3 b / 2 



<>p j t-> ' 

where C (3 / 2) = ^ ~ 19. 

The dependence of Eq.(^) is shown in Fig. 4. Therefore, any planet 
(even consisting from pure hydrogen) should have a mass less than 10 31 g (if 
its density is approximately equal to lg/cm 3 ). 

In Fig. 4 the masses of the planets of the Solar system are marked. The 
mass of the Jupiter is 1.9 ■ 10 31 g. It is close to the specified limit. For the 
Jupiter Eq.(p8[) gives (5 ~ 2. It is according to the data that the large planets 
have the deuterium-helium composition. For other planets the mantle is not 
small in comparison with their sizes. For this reason, Eq. (|28|) can give an 
excessive estimation for other planets. 
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Figure 4: The dependence of the core mass of planets on (3 (Eq.(|28|)) at 
7 = lg/cm 3 . On the ordinate, the logarithm of mass (over lg) is plotted. 
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6.2 Relativistic electron-nuclear plasma. 

When the pressure increases, the substances are transmuted into relativistic 
electron- nuclear plasma (the polytropy k= 3). Its state equation is [|7| 



PW = : 4/s ' (29) 



( 37T 2 ) 1/3 ftc 7 4 / 3 

4mp /3 /? 4 /3 

If this plasma is originated by the GIEP effect, then the steady-state 
value of mass of a star consisting of it, according to Eqs. (^2"lj2"9D is 



M (3) = C (3) ■ AT ■ ^ (30) 



where the dimensionless constants are 



A * = fe) = 154 ■ 1038 (31 > 



and G( 3) = (1.5 5 tt) 1/2 ~ 4.88. 

Because of the electric neutrality, one proton should be related to elec- 
tron of the Fermi gas of plasma. The existence of one neutron per proton is 
characteristic for a substance consisting of light nuclei. The quantity of neu- 
trons grows approximately to 1.8 per proton for the heavy nuclei substance. 
Therefore, it is necessary to expect that inside stars 2 < (3 < 2.8. 

The masses of stars can be measured with a considerable accuracy, if 
these stars compose a binary system. There are almost 200 double stars 



which masses are known with the required accuracy |T7| . Among these stars 
there are giants, dwarfs, and stars of the main sequence. Their average 
masses are described by the equality 

(M star ) = (1.36 ±0.05) M , (32) 
where M is the mass of the Sun. 

The center of this distribution (Fig.5) corresponds to Eq.(^) at (3 ~ 2.6. 

6.3 Ultrarelativistic electron-nuclear plasma. 

Further increase in pressure transmutes substances into ultrarelativistic plasma. 
Then nuclear reactions of capture of electrons by nuclei become favorable and 
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Figure 5: Mass distributions of stars and pulsars from the binary systems 
[|TH|. The curve shows (3 (Eq.([HJ)). 
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the neutronization of matter takes place. Equilibrium pressure of ultrarela- 
tivistic plasma does not depend on its density. It is formally characterized 
by the polytropy k=-l and its state equation is [0] 

A 4 , , 

p (-V = 19 i77 x3 - ( 33 ) 
lz7r 2 {ncj 

Here A is the difference between the energy of the initial nucleus and the 
energy of the daughter nucleus. 

The equilibrium mass of a star, consisting of ultrarelativistic plasma, 
according to Eqs.(p^),Eq.(|33"D is 



M ( _i) = Cr-D ( ^— ) , (34) 

1 } 1 ' \(hc) 9/ 2 G 3 / 2 7 2 / 



whereC ( „ 1) = i i,(| 7 ) 1/2 ^6-10- 3 . 

According to the astrophysical data a neutronization of matter takes place 
at the density 7 lO 7 ^. Thus, Eq.(|34|) gives 

Af ( _i) « 8 ■ 10 32 # « O.4M . (35) 

Certainly this result is the rough estimation on the order of magnitude 
only, but it is in a satisfactory agreement with measurements of the as- 
tronomers related to masses of white dwarfs from double systems. 

6.4 Nonrelativistic neutron matter. 

At higher pressure, the substance is transmuted into a nonrelativistic neu- 
tron matter with a small impurity of protons and electrons J7|. The state 
equation of the nonrelativistic neutron matter will coincide with Eq. ( p2|) with 
a replacement of m e with m p 

P|3/2, ~ SO 3 ' 1 ' 

Together with Eq. (p2|) , it gives the equilibrium mass of the nonrelativistic 
neutron star 
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«&> = cm(ty i^rr (37) 



As the density 7 ~ 4 • 10 13 g/cm 3 and (3 = 2.6 

M$ 2) ~ O.O5M . (38) 
The astronomers have not found such neutron stars. 

6.5 Relativistic neutron matter. 

With further increase in pressure, the neutron Fermi gas becomes a relativis- 
tic one. Its state equation completely coincides with the state equation of 
the relativistic Fermi gas of electrons Eq. (p2|) 

As well as the masses of relativistic stars, the masses of relativistic pulsars 
do not depend on their density and can be directly expressed through world 
constants 



M$ = C {3) ■ AT ■ % (40) 



J 2 

(at /3 — 1 for the pure neutron Fermi gas). 

As it is specified in [0, at this density of matter it is necessary to take 
into account nuclear forces and the presence inside nuclear matter except 
neutrons also of protons, 7r~mesons, and electrons. It can be made using (3 
as a parameter of the correction. 

The mass of the neutron star can be measured with a considerable ac- 
curacy if it enters into a binary system. The astronomers have found 16 



radio-pulsars |L8[ and 7 X-pulsars |16| in binary systems. They all are lo- 



cated in a very narrow mass interval (Fig. 5) 

(M pulsar ) = (1.38±O.O3)M . (41) 
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The center of this distribution corresponds to Eq.(f40D with the correction 
parameter (3 ~ 2.6 just as for relativistic stars. Thus, we come to the conclu- 
sion that Eg. ([TOD on the order of magnitude correctly describes the results 
of astronomical observations. 

7 Conclusion. 

It is expedient to underline the basic obtained results in summary. 

1. The developed theory defines a concept of the steady-state values of 
masses of celestial bodies related to their equations of state and gives the 
possibility to calculate these values. 

2. It gives the new way for the determination of the substance density dis- 
tribution inside celestial bodies. According to early models, it was supposed 
that density of a substance inside celestial bodies grows more or less mono- 
tonically with depth and at the centre of a star, the density has the greatest 
value and even a black hole may exist there. According to the developed 
theory, the density of a substance inside a star is constant. 

3. It is interesting to emphasize that the "biography" of such a star 
appears much poorer than in the Chandrasecar model. 

There cannot exist a black hole inside a star, and it should not collapse 
with a temperature decrease. All the considered effects are based on an equi- 
librium of the Fermi system. Temperature does not influence the parameters 
of relativistic plasma. Therefore, a star with a mass close to the steady-state 
value (Eq. (|30D ) is in a stable equilibrium not depending on temperature. The 
existing stars should retain the stability at any (even at zero) temperature. 
Therefore, a collapse of the already existing stars apparently is impossible. 
The instability of a star can arise with burning out of light nuclei - deuterium 
and helium - and with a related increasing of (3. This growth leads to the 
reduction of a steady-state value of mass (Eq.([H])) and, probably, to the 
distraction of the stars wiht masses more than the steady-state value. 

4. The developed theory determines the simple and essential mechanism 
of generation of the magnetic field by celestial bodies. All early models tried 
to solve the basic problem - to calculate the magnetic field of a celestial body. 
Such a statement of the basic problem of planetary and stellar magnetism 
is unacceptable at present. Space flights and a development of astronomy 
discovered a remarkable and earlier unknown fact: the magnetic moments 
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of all celestial bodies are proportional to their angular momenta and the 
proportionality coefficient is determined by the ratio of world constants. The 
explanation of this phenomenon is the basic problem of planetary and stars 
magnetism nowadays. Early models cannot explain this phenomenon. The 
developed theory used for this explanation a standard mechanism. 

5. It is possible to consider that now the predicted steady-state values of 
masses of celestial bodies and the predicted values of their magnetic moments 
are in a satisfactory agreement with the data of observations. But it is 
tempting to obtain these data to closer limit of accuracy. Two arrows in the 
upper part of Fig. 5 mark masses of stars consisting of extremely light and 
heavy nuclei. These values are obtained from Eq . (|30D without the use of any 
fitting parameters. In agreement wiht the developed theory, if stars have a 
"usual" chemical composition, there must be no stars outside of this interval 
(or these stars should be unstable). The histogram on Fig.5 is somewhat 
wider. It is interesting to understand, whether there is a principal deviation 
from the developed theory or it is a result of measuring errors. First of all, it 
requires a more careful and precise measurement of masses of binary stars. 
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